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Abstract 

A fluid of Skyrmions coupled to the dilaton field and the w me- 
son field is considered. A mean field theory is developed in which 
the dilaton and the lo field acquire a mean value determined by the 
Skyrmions. The influence of the background fields on the Skyrmion 
profile is determined and consequently the scaling properties of the 
Skyrmion follow. The model obeys chiral invariance and scale invari- 
ance broken only by the dilaton trace anomaly potential. The dilaton 
plays the role of the scalar field as in the a + lo model. The dilaton 
potential is augmented by terms that do not spoil the trace anomaly 
while allowing a fit to nuclear matter properties. The phase diagram 
of the model shows unique features, like a lack of solutions for certain 
densities and temperatures, signalling the appearance of a new phase 
that can not be described in terms of Skyrme baryons. 

PACS 12.39DC, 21.10Dr, 21.65.+f, 24.10.Pa 
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1 Introduction: Application of Skyrmions to 
the many body problem 

Around 35 years ago Skyrme |l|] proposed to describe baryons as topological 
solitons in a meson field theory. More recently, studies of the large Nc approx- 
imation to QCD PI suggested that indeed the meson degrees of freedom are 
dominant and baryons might arise as solitons @, H- In the last decade there 
has been a large body of research on these solitons, or Skyrmions @, ||, H, 0]. 
The research has mainly focused on single baryon and few baryon systems. 

Skyrme himself was interested in applying his model to systems of many 
nucleons aimed at the description of nuclei. After gaining support from QCD 
it appears that a model of nuclei based on Skyrmions would be better related 
to more fundamental theories of matter. This is especially appealing due to 
the difficulty of solving QCD explicitly. 

The Skyrme model has had moderate success in dealing with the nucleon- 
nucleon interaction, including the attractive isoscalar central potential for 
which several mechanisms were proposed. It appears then natural to in- 
vestigate more formal issues involving Skyrmions, such as their behavior in 
nuclear matter. 

In the present work a mean field model of dilatons and u mesons coupled 
to Skyrmions is investigated. The model is inspired in the classical a + 
uj Walecka approach [B. The original a + u model has been extended in 



order to include scale invariance by means of the dilaton field in refs. [pi) 
pj| . In these works two scalar fields are used, the chiral partner of the pion 
represented by the a field and the dilaton. The models try to reproduce 
correctly both the bulk properties of nuclei and their excitation spectra. 
Chiral symmetry is introduced in its nonlinear realization in the former and 
in its linear one in the latter, while the nucleons are considered as pointlike 
Dirac objects. A different attempt to include chiral symmetry in the mean 



field theory of Walecka was done by Lynn IT^. The virtue of the Skyrme 
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model is that it enables to include both chiral symmetry and baryons in a 
consistent way, the baryons being solitons of finite extent in space. These 
aspects make the model a good tool for the investigation of baryonic nuclear 
matter, but at the same time the treatment becomes more involved. We will 
see, however that in a dilute fiuid approximation it is possible to treat the 
Skyrmions as essentially free particles interacting with a medium that carries 
the information of density and temperature. This medium acts through the 
mean fields of the dilaton and the uj that are determined by the Skyrmions 
themselves. Mishustin [|13[ has treated a single Skyrmion immersed in a bath 
of Dirac particles. Our approach is to generate the bath from the Skyrmions 
themselves in a self-consistent manner. Section 2 deals with the dilute fiuid 
approximation. Section 3 describes a specific Skyrmion fiuid. Section 4 treats 
the choice of dilaton potential. Section 5 shows the results of the mean field 
approximation for finite temperature and density. 



2 Dilute Skyrmion fluid 



Consider a field theory lagrangian of Skyrmions, the dilaton a and the uj 
meson 0, [l^ 



Here 



2 dilaton 



+ C2 + Ci-V, 



interaction 



•2a 



}-Tld,ad^a-^iT{L,L^ 



- V{a) + 



B[l + e^'^(4a - 1)] - \{d,u. - d.u.f + \e'^ml (1) 



= U^d^U, (2) 

where U{r, t) is the chiral field, is the pion decay constant and e the 
Skyrme parameter. 

The trace of the energy-momentum tensor in QCD is given by 
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TH = d,D^ = -^G^C^^^ ^ (3) 

where D^{= T^^x^) is the dilatation current, as is the QCD couphng con- 
stant, is the gluon field, and tp is an order-parameter field — the dilaton — 
which represents the scalar glueball formed from the contraction of the two 
gluon fields; is taken in the fourth power in anticipation of its use as a 
scalar field of dimension 1. In the lagrangian of eq. (|1]) ^/^ ~ Tqc'^. The trace 
of the energy momentum tensor of eq. ([^) -generated by V^- is constructed 
in accordance with eq. (|). 

The lagrangian incorporates scale invariance broken only by the anomaly. 
All other terms are scale invariant, including the uo mass term that carries 
appropriate factors of . 

In the present work we focus on the SU(2) version of the Skyrmion. The 
lagrangian possesses chiral as well as isospin symmetry. We have omitted a 
pion mass term that manifestly breaks chiral symmetry because it is relatively 
unimportant. We have also omitted p meson coupling as it is our intention 
to treat symmetric nuclear matter only. The two parameters Fq and B are 
related to the glueball — or scalar field — mass through = IG-B/Fq. An 
accepted value of the mass is around 1.5 GeV, although the mixing of the 
glueball with vrvr-exchange may bring it down to the below GeV range. (A 
critique to the use of the dilaton in nuclear matter applications can be found 



in ref. ITS 



Analogously to the mean field theory of pointlike baryons, we consider an 
ensemble of essentially free Skyrmions. Each Skyrmion will be accompanied 
by its own dilaton and uj fields. Although we deal with free Skyrmions, the 
average properties of the ensemble are still included in distribution functions 
that depend on the density and temperature. 

For the N Skyrmions we use the product ansatz 



4 



UB=N{r, Ri, R2, ■ ■ ■ , Riv) = U{r - Ri) U{r - R2) ■ ■ ■ U{r - R^v), (4) 
whereas for the scalar fields we use an additive ansatz 



(rB=N = 0"l + 0-2H hfTAT, 



aQ + 6ai + 6a2 + h SaN, 



(5) 



UJb=N = LUi + UJ2 + ■ ■ ■ + ^N, 



UJq + 6uJi + 6uJ2 + ■ ■ ■ + Sujn 



(6) 



Where cro,ciJo are the mean field constant values of the dilaton and the 
uj and 5cr, 5uj represent the fiuctuations. The fields a, 5cr, cj, 5uj depend on 
the same arguments as the Skyrmion they are attached to. The ansatze 
above insure total baryon number = N and allow an easy separation of the 
mean field degrees of freedom from the microscopic excitations inside each 
Skyrmion. In thermal equilibrium, the mean field fields will depend on the 
temperature T and the chemical potential /i. Self-consistency then requires 
that the values of do and ojq have to be determined by the properties of the 
ensemble. For a certain distribution function f(yU, T) the temperature, we 
then demand the phase space thermal average to be 



with a similar equation for uo. The mean field average value of the fiuctua- 
tions vanishes by definition. In eq. the integration is over the collective 
coordinate coordinates Ri of the Skyrmions and their conjugated momenta. 



< 0" > = < ctq + 6ai + ■ ■ • + da^ > 




(7) 
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In the diluted fluid approximation, Skyrmion interactions are neglected and 
consequently there are no potential terms depending on the coordinates Rj 
of each baryon. There remain only kinetic energy terms for each Skyrmion 
individually. Standard quantization of these terms will then determine the 
corresponding wave functions to be plane waves, as for the conventional Fermi 
gas model. 

A crucial point in the mean field theory of skyrmions is the topological 
baryon density 



B^" = ^tr 

247r2 



(8) 



where e'^"^''' is the totally antisymmetric tensor density, and it is easily seen 
that 

d^B^ = 0. (9) 
The corresponding baryon number becomes 



Bo = ^tr 

" 247r2 



U^d.u) {U^d.u) (u^dkU)] , (10) 



Using the product ansatz of eq. (|^) in the baryon number above allows 
us to calculate the mean field average by integrating over the Skyrmions' 
coordinates and momenta. Consistently with the approximation of a fluid of 
free Skyrmions we demand that there is no overlap between the Skyrmions 
profile functions, and further assume that the profile drops to zero at a dis- 
tance smaller than the inter-Skjnrmion separation. This condition is needed 
in order to insure integer baryon number within each Skyrmion cell, because 
the single baryon density depends on the values of the profile function at the 
particles' location and at infinity. In the present case, infinity is replaced by 
a finite distance. Typically, this distance is of the order of 0.8 fm, whereas 
the interparticle separation in normal nuclear matter is approximately 2 fm. 
It is then not very unrealistic to assume that most of the topological baryon 
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density is contained within a volume smaller than the average volume avail- 
able to each Skyrmion. Note that, even if there is some overlap between the 
Skyrmions, the choice of product ansatz insures total baryon number = A^. 
To this approximation we obtain 

Bq = 61 + 62 H \-bN, 

<Bo> = N/V (11) 

where bi is the baryon density of the i^^ Skyrmion, V is the volume of the 
fluid, and the average is defined in eq. (^. The integration over the collective 
coordinates of the Skyrmions is mathematically equivalent to an integration 
over the coordinates r in eq. (^, this is the reason why the averaging above 
essentially counts the number of Skyrmions. After the Fermi averaging, the 
localized baryon density becomes uniformly smeared out. There is no residue 
of the individual Skyrmion densities. We will see below that the influence 
of the mean fields on the Skyrmions operates solely through the dilaton 
field. Although the baryon density of each Skyrmion will vary accordingly, 
the integrated density is scale invariant and will still contribute the value 
B = 1 for each Skyrmion. The dilute fluid approximation is therefore closely 
related to the mean field Dirac model with two major differences: 1) there 
is a different dynamics for the baryons dictated by the lagrangian of eq. (|^) 
and, 2) there is a clear way to uncover the baryon response to the medium. 

3 Mean field Skyrmion fluid 

The many body Skyrmion fluid in the dilute approximation is built from the 
product ansatz of eq. (^) for the Skyrmions and the additive ansatz for the 
mesons of eq. (|^). We therefore first focus on the single baryon case. Each 
soliton is constructed as a static solution of the equations of motion. The 
collective coordinate energy of the Skyrmion translation is later introduced 
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in an adiabatic approximation. The single baryon ansatze are then 



f/(r) = exp[iT ■ rF{r)], 
u^{r) = (u;(r), 0,0,0) 



(12) 



Substituting eq. ([T^ ) into eq. ([^) we find the static mass of the Skyrmion 

M 



M{r) 



A-n / r^drM{r) 

sin^F 



+ 2- 



+ 



1 sin^F 

2g2 y,2 



sin^F 



+ 2F 
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+ V{a) 



I 2cj 1 7^2 ^/2 1 , ,/2 I 



Qv uj F' sin F 1 2 2 2<t /1Q^ 
r— m , UJ e 13 



where primes denote derivatives with respect to r. The Euler-Lagrange equa- 
tions for the static profile and the mesons become 



^2. Ssin^FN^,, 2^^,/ / r. 4sin2FF'2 e2<^sin2F 

e^"" + — - — F" + 2e^'"F'(l/r + (j') 



^2 / ^2 ^2 



4sin2 F sin 2F ^ 2^y u;' sin^ ^ _ q 



.^2,^2 .^2 



Tl e^" o" + 2a'2 + 2 F'^ + 



(icT 

" , 2 2<T , gyF'sm^F 

0;+ - m^uj e + — -^-^ — = 14) 



where f = eF^rr and again primes denote derivatives with respect to r. 

Using eq. (|^) and the fact that the fluctuation of a and lo vanishes in the 
mean field state we readily find that except for small contributions coming 
from the a potential the equations above are modified in the presence of a 
nonvanishing cuq and ctq by the simple scaling law 

r ^ e~'"Or , w ^ e'^Ocu (15) 
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It is then unnecessary to solve the equations of motion for the single Skyrmion 
with the meson fluctuations included. It suffices to solve them for a free 
Skyrmion and then rescale the uo field and the radial distance. Moreover, 
for dilaton masses of the order of 1 GeV and more, the fluctuation 5(j is 
altogether negligible. Therefore the scaling law becomes exact. Analogously 
it is easy to show that the static mass of eq. (|^) scales as 

M = Mqc"" (16) 

where Mq is the mass for o" = 0. 

Before proceeding to the many body case, we consider the Lorentz boost- 
ing of the static Skyrmions. We perform a boost on the collective coordinate 
R{t) of each Skyrmion. For the sake of simplicity consider a Lorentz boost 
along the x axis with velocity parameter v. 



a; — > X 



X - R{t) 



y = y 

z = z 

F{r) F{r) 

(T(r) ^ a(f) 



with uJ^B^ = u;o(f)i?o(f). Introducing the above transformation in eq. (j^) 
and calculating the Hamiltonian, we find the energy of a Skyrmion in motion 
to be 
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where 



y + M2, p 



Mv 



M is the static mass of eq. (|T3p for nonvanishing ujq, ctq and 



E2 
E^ 

Uint 



477 



sin^F 



1 r^rfr e^V'^ 

An J r'^dr K 

An / r^c/r e'^" ^ — 



TT 



J dr uj F' sin^ F 



(19) 



The energy Ep of eq. (0) enters the single particle distribution functions 



rir 



^exp 


"ep - jj,' 
[ kT \ 


-) 


^exp 


[ kT \ 





-1 



-1 



(20) 



where ep = Ep + gv ujq, tp = Ep- qvujq. 

We can now write down the energy of N Skyrmions per unit volume in 
the mean field approximation for symmetric nuclear matter. 

——3 Ep{np + Up) + \4(fTo) - - e^'^om^^o + 9v i^o Pv (21) 
(2 tt) 2 



Ev = A 
where 



Pv 



d^p 



3 K^p ^p) 



(22) 



At T = we have Up = Q{pf — p),np = yielding the equations of motion 
for the mean fields 
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dEv , f d^p dEp dVa 2 2^0 2 



At finite temperature we have to maximize the pressure instead of minimizing 
the energy. The contribution to the pressure of a single Skyrmion is given by 

= [E. + E, - B„)^^5^ + ^^-^ --((/,- (24) 
and the pressure of the ensemble per unit volume becomes 

Pv^^j ^ PAn, + Up) - V,{a^) + ^e^^»m>o' (25) 

As mentioned above, due to the large dilaton mass one can neglect the 
fluctuation of the dilaton field inside the Skyrmion. Numerical evaluation 
of 5a for a dilaton mass above 1 GeV indeed gives a negligible dilaton field 
5a ~ 0. Moreover, the simplest approximation to the mean field demand 
of a vanishing expectation value of the fiuctuation of the uj field is 5oo = 0. 
This first approximation to the problem permits a large simplification of 
the single particle's energy and pressure. Using the virial theorem for the 
soliton profile appropriate for the case of vanishing dilaton and cu fluctuation, 
E2 — E4 — |M , we obtain 

M = nFy^^ I r^dr [f'^ + ) = e'^Wo 

= 4 ^^'^ 

To this approximation the mean field Skyrmion fiuid becomes identical to 
the Dirac mean field approach, with the bonus of knowing how to calculate 
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the reaction of the single Skyrmion to the bath using the scahng properties 
determined by the dilaton. As it is our goal to investigate the most important 
contribution to the Skyrmion properties due to the bath, we will proceed 
along this line and neglect the dilaton and the uj meson fluctuations inside the 
Skyrmion. The results we will obtain will then be applicable to a conventional 
mean field theory of Dirac pointlike particles coupled to the dilaton and the 
u meson. The key new ingredient in such a treatment would then be the use 
of a single scalar field in order to fit the nuclear matter properties. 



4 The dilaton potential 

In eq. (|l|) we introduced the conventional glue potential that reflects the 
trace anomaly |T6|, where the 'bag constant' is 5 ~ (240 MeV)^. The po- 
tential can be supplemented by other terms that do not affect the anomaly. 
Modifications of the potential have been proposed in order to include chiral 
symmetry breaking corrections and quark contributions to the trace anomaly 
[ P^l , whereas Furnstahl et al. |]TU| opted for an approximation to the potential 



that takes into account anomalous dimensions acquired by the scalar fields 
upon renormalization. In any event there is a need to modify the potential in 
order to fit the properties of nuclear matter. The difficulty arises mostly in 
fitting the low value of the nuclear compressibility modulus n ~ 270 MeV. In 
the abovementioned works two scalar fields come into play, the dilaton and 
a scalar partner of the pion, the a meson. In the present approach there is 
no room for the latter, because we are working in the nonlinear realization of 
chiral symmetry due to the use of the Skyrme lagrangian. We therefore allow 
for modifications of the dilaton potential that do not affect the anomaly, but 
influence the predictions for nuclear matter. Terms of the form e"'^ — 1 , 



generically referred to as no-log terms |TT[, can be added to the potential 
leaving the vacuum expectation value of the a to be that dictated by the 
anomaly potential, chosen here to be cr = by a simple shift of the original 
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dilaton field. The actual choice of the functional form of the new terms in 
the potential turns out to be relatively unimportant. (There exists also the 
possibility of modifying the uo potential as in ref. |]1T|], introducing higher 
order terms in the uo field). 

Let us consider the constraints on demanded by nuclear matter phe- 
nomenology. The potential has to be such that at the saturation density of 
nuclear matter po = -154 baryons/fm^: a) the binding energy per nucleon is 
16 MeV ; b) the binding energy is maximal ; c) the compressibility is of the 
order of 270 MeV; d) the dilaton and uj fields obey the mean field equations. 
These conditions translate into the following equations 



gy Po e 
ml 



-2<70 



UJQ = 



dV„ 



- ulmy-' + Q 







(27) 



where 






Po(Mo - 16 MeV) 




(28) 




(29) 
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where 



d'Ey ^ pI ^ gl 



In order to fulfill the above constraints we introduced an additional piece 
in the dilaton potential with three new free parameters that proved to be 
quite successful in reproducing the nuclear matter phenomenology 



V^ = B[l + e''^iAa-l)] + B 



ai(e-"-l)+a2(e"-l) + a3(e2--l)+a4(e=^"-l 

(31) 

where B is fixed and the anomaly condition requires 

^ = 

da 

(32) 

at 0" = implying Oi = a2 + 2a3 + 3a4. 

We have chosen the term multiplied by ai to have a negative power of a 
in order to avoid the introduction of a second minimum in the potential for 
a < 0. The only sensible minimum then remains the one at a = 0. There is 
a need to introduce three new parameters in order to reproduce correctly the 
nuclear matter binding energy at the right density, and the compressibility 
factor. 

The potentials for the fluctuations of the a and u fields are then deter- 
mined by the averages of eq. (|^) to be to lowest order a massterm interaction 



1 



2 
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V{6a) = B-Sa^ 16e^^° (1 + 4ao) + aie"''" + 026'"° + 2a3e^'"' + a^-e 



3ao 



2 



14 



V{6uj) = ]^^''^'mlbu? (33) 

The effective mass of the dilaton is not fixed, because it depends on the 
parameter Fq of eq. (|^), that has to be prescribed separately. The mass will 
enter only in finite nuclei calculations. 



Results and discussion 



The conditions on the dilaton potential of eqs. , |29|) were implemented 

by choosing a suitable effective mass at the nuclear matter saturation density 
and a certain compressibility. We took M* = 694 MeV corresponding to 
ao = —0.3 and k = 300 MeV. The value of k we used is very close to the 



measured one [0 - although much lower values are suggested by the analysis 
of Blaizot [|I9]- ,whereas for the effective mass we took a conservative figure, 
although many authors tend to select a lower value of around 600 MeV. We 
have found that there are no qualitative differences in the predictions of the 
model when the latter mass is used. In all the fits we kept the bag constant 
at B = (240 MeV)"^, although we were able to fit the data wit a large range 
of values of this parameter. 

With the above choice we find : ai = -5.53, a2 = -54.74, as = 48.63, 04 
= 16.01 and, gy = 7.29. Although it may seem that the dilaton potential 
coefficients are a bit large, one has to bear in mind that they multiply ex- 
pressions very close to zero for all the attainable values of a. The u meson 
coupling constant turns out to be smaller than expected but still within the 
range quoted in the literature. Other terms in the u meson potential might 
be needed in order to allow higher values of gv- 

Figure 1 shows the binding energy per nucleon as a function of the density 
both for the normal and abnormal solutions. The normal branch possesses 
all the desired requirements for nuclear matter, therefore at T = it rep- 
resents the physically stable phase. The existence of an abnormal branch 
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is well known. It arises from the functional dependence of the scalar field 
potential. This branch is characterized by large negative values of the scalar 
field. In some models the normal and abnormal branches cross each other 
at some critical density, suggesting a transition to a dense, chiral symmetry 
restored phase, accompanied by a dramatic decrease in the nucleon effective 
mass. Chiral symmetry is built-in in the Skyrmion fluid, hence, the effec- 
tive nucleon mass is not an indicator of chiral symmetry restoration. The 
abnormal branch is here totally spurious, and might be possibly reached in 
a metastable situation. 

At finite temperature the criterion of maximal pressure for a fixed chem- 
ical potential determines the physical phase. Again, the normal branch 
prevails. 

Figure 2 shows the nucleon effective mass as a function of density and 
temperature. (The apparent cusp in figure 2 is due to the rudimentary plot- 
ting procedure). 

In the Walecka model -and many others similar to it- the nucleon mass 
decreases as a function of density This is not the case in the Skyrmion 
fluid model. The minimum effective mass arises here at a density of around 
1.7 po- Contrarily, the usual parametrizations of the scalar field potential in 
Walecka models, do not show this limitation. The reason for this difference 
can be traced back to the dynamics dictated by the dilaton, especially to the 
modified trace anomaly potential. 

The scalar field potential of the Walecka type models is built in terms 
of powers of the scalar field, up to fourth order generally, to satisfy renor- 
malizability requirements. (Although this condition is usually hard to meet 
due to the need for a negative coefficient for the fourth order term in most 
successful fits.) However, the dilaton potential is made to obey the condition 
demanded by the trace anomaly. Only no-log terms that do not spoil the 
anomaly are allowed. This strong demand makes the potential qualitatively 
different. The solutions to the nuclear matter equations with this potential 
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for the normal branch, exist only in a narrow range of values of a near a = 0. 

It is now possible to understand why the effective mass increases beyond 
a certain density: The dilaton attractive contribution to the mass is limited, 
while the uo meson repulsion grows in direct proportion to the nuclear matter 
density, eq. (p7|). The solution of the nuclear matter equations then tend 
to push the dilaton towards positive values in order to fulfill the scaling 
properties of the model, eq. (|15|) . 

The Skyrmion properties are essentially determined by the scaling char- 
acteristics of the model. The mass of figure 2 scales as e""" whereas the root 
mean square radius scales as e~"° . The Skyrmion swells as a function of den- 
sity up to densities below 1.7po- This is due to the attractive effect induced 
by the dilaton. Above a certain density the Skyrmion starts to shrink, due 
to the abovementioned repulsive effect in a manner independent of temper- 
ature. This is reminiscent to the behavior of the Skyrmion interaction at 
short range (high density) that is known to be repulsive and strong. Other 
properties, like qa-, magnetic moments, etc, follow directly from the scaling 
of eq. (|1|). 

Figures 3-5 depict the pressure, chemical potential and entropy as a func- 
tion of density and temperature. Although it may appear that there could 
be a phase transition between T = 125 MeV and T = 175 MeV, calcula- 
tions performed with a finer grid than the one shown in the graphs, does 
not support such a scenario. The transition between those temperatures is 
smooth. 

At a temperature of 25 MeV the graphs of the pressure and chemical 
potential show a rough behavior due to the proximity of this temperature to 
the liquid-gas phase transition that appears in a broad class of Walecka type 
models 0. 

The entropies of figure 5 were calculated with three different methods. 
The agreement between the numbers obtained determine the accuracy of the 
minimization procedure to be better than 0.5 %. 
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As shown in figure 6, above T = 190 MeV there cease to exist solu- 
tions around the the normal density of nuclear matter. Above T = 220 MeV 
solutions cease to exist for densities below 4po, perhaps signalling the appear- 
ance of a different phase that we can not describe in terms of Skyrmions and 
mesons. Other degrees of freedom must come into play, perhaps quarks and 
gluons. This behavior is different as compared to the results of conventional 
models [|l7l for which there are solutions at all densities and temperatures 
with a phase transition to a baryonic plasma phase above a certain temper- 
ature typically of the order of 200 MeV. The reason for this difference can 
be traced back to the special features of the dilaton coupling and potential. 

In the present work we have paved the way for a treatment of a Skyrmion 
fluid based on the principles of scale invar iance and chiral invariance. In a 
later work we will test the approach and especially the potential obtained 
in finite nuclei calculations in the Thomas-Fermi approach as well the cor- 
rections due to the rotational degrees of freedom of the Skyrmion. In a 
finite nucleus calculation new parameters will enter, like the dilaton mass. 
It is then a bit risky to venture any predictions for such a system. A much 
more complicated and ambitious task would be to tackle the introduction of 
Skyrmion interactions. 
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Figure Captions 

Fig. 1: Binding energy per nucleon in MeV as a function of p/po- Normal 
branch, dashed hne and abnormal branch, full line. 

Fig. 2: Effective Skyrmion mass M* in MeV as a function of p/po- Full line, 
T = 25 MeV, dashed hne T = 75 MeV, dash-dot line, T = 125 MeV 
and dotted hne, T = 175 MeV. 

Fig. 3: Pressure in MeV/fm^ as a function of p/po- Line definitions as in Fig. 
(2). 

Fig. 4: Chemical potential in MeV as a function of p/po- Line definitions as in 
Fig. (2) 

Fig. 5: Entropy per baryon as a function of p/po- Line definitions as in Fig. 
(2) 

Fig. 6: Pressure in MeV/fm^ as a function of p/po- Dotted line, T — 190 MeV, 
fuU line, T = 200 MeV and dashed hne, T = 220 MeV. 
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